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ABBREVIATIONS AND SYMBOLS 



ROMAN SYMBOLS 
A Area 



'D 



K-e model empirical constants (Table 1) 



Cp Specific heat at constant pressure 

E 9.0 

h Enthalpy 

h Stagnation enthalpy 

I Turbulence intensity 

K,k Turbulence kinetic energy 

M Average molecular weight 

P,p Pressure 

r Radial distance 

R Gas constant 

R Universal gas constant 

5 Source terms 

T Temperature 

u Axial velocity 

V Radial velocity 

X Axial distance 

yp'*’ Dimensionless distance from solid boundary 
GREEK SYMBOLS 

r Effective transport coefficient 

6 Incremental distance from wall 
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e Turbulence dissipation rate 

K von Kantian constant 

y Viscocity 

p Density 

a Prandtl or Schmidt Number 
T Shear stress 

<p Any variable 

SUBSCRIPTS 

eff Effective 

lam Laminar 

p Near wall node 

t Turbulent 

w wall 
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I. 



INTRODUCTION 



During the past few years, there have been many advance- 
ments in numerical techniques for predicting the behavior of 
fluid flows. For example, several computer models have been 
developed by Gosman, Spalding and others [1,2,3] which use 
the mass, momentum and energy conservation equations reduced 
to finite difference, nonlinear algebraic form. The develop- 
ment of reliable computer programs of this type greatly bene- 
fits engineering analysis in such widely varying fields as 
meteorology, aerodynamics and gasdynamics. 

The earlier two-dimensional computer codes were based on 
vorticity (co) and stream function (iff) [1,2,5]. This form of 
the governing equations eliminates pressure and velocity from 
immediate consideration. Pressure is normally calculated 
only after a converged solution is obtained. This technique 
has several inherent disadvantages: 

1. It results in large errors in the predicted pressure 
distributions in all but quiescent flow regions due to the 
higher order dependence of the pressure gradient on stream 
function [6], 

2. It is usually restricted to constant density flows 
or to flows in which density varies only with temperature 
[3,6] . 

3. The boundary conditions are difficult to specify 

[3,5] 

4. Considerable difficulty has been experienced in 
arriving at converged solutions, especially for non-uniformly 
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spaced grids and high flow rates [2, 5, 6, 7]. 

5. The model is not easily extended to three- 

dimensional flows [3]. 

To overcome these difficulties, emphasis has been placed on 
developing computer codes based on velocity and pressure, the 
primitive variables. 

A major problem with any new computer model is model 
validation. The difficulties of collecting accurate enpirical 
data are multiplied when investigating three-dimensional and/ 
or reacting flows. In addition, many variables within these 
flows are not readily measurable (turbulence intensities, etc.). 

An effort to utlize elliptic conputer models which can 
handle turbulent, reacting, variable density flows at high 
subsonic and sonic velocities has been underway at the Naval 
Postgraduate School for several years. Two specific areas 
which have been investigated are flows in a turbojet test cell 
and in the combustion environment of a solid fuel ramjet. 

It is important to have the capability to test high per- 
formance jet engines throughout their operating envelope 
under conditions which approximate installed conditions. This 
is often accomplished in blockhouse type installations called 
turbojet test cells (TJTC) . The typical test cell incorporates 
an inlet, a horizontal test section and a vertical exhaust 
stack. The engine to be tested is normally mounted near the 
center of the cell to allow the development of a nearly uniform 
engine inlet velocity profile. The engine exhausts into an 
augmentor ti±ie which entrains additional air for exhaust gas 
cooling and dilution. The quantity of this secondary air is 
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crucial to proper engine testing and test cell performance. 

Testing today’s high power and high mass flow engines in 
these installations produces a myriad of noise and air pollu- 
tion problems. Cell modifications must often be made to mini- 
mize these problems. This fact coupled with the future need 
for larger, more expensive test cells to replace obsolete cells 
and to accommodate new generations of high technology engines, 
makes the development of reliable modeling methods imperative. 
The frequently used one-dimensional models are not adequate 
for predicting the details of the complicated flows within a 
turbojet test cell and, therefore, the cells often do not per- 
form to their designed limits. An accurate flow model would 
provide a needed design tool which could help prevent costly 
design errors and improve operating efficiency. 

A two-dimensional rp-co computer code was used to analyze 
the flows in a full scale and a subscale turbojet test cell 
at the Naval Postgraduate School [2,6]. Experimental data 
from the subscale test cell have been compared with computa- 
tions made with this computer model [6]. A primitive variable 
(u-v-p) computer model could improve this capability by ex- 
tending it to specific geometries and flow rates that the 
i])— oj model was incapable of predicting. In addition, a u-v-p 
model would more readily allow variable density flows to be 
analyzed and should more accurately predict augmentor pressure 
distribution. An application of a primitive variable model 
to jet pumps has been reported by Croft and Lilley [4]. 
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The primary objective of this investigation was to adapt 
a primitive variable model to the turbojet test cell geometry 
and flow conditions and to validate that model with data from 
a subscale turbojet test cell. Possible utilization of the 
model for predicting optimum placement of sound suppressors/ 
jet break-up devices within the augmentor tube was also inves- 
tigated. 
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II . MODEL OVERVIEW 



A. INTRODUCTION 

The computer model used in this study was adapted from 
the CHAMPION 2/E/FIX computer program developed by Pun and 
Spalding [8]. CHAMPION is a TWO- dimensional Elliptic, Fix ed 
grid conputer program which provides a solution of the conser- 
vation equations for recirculating flows in finite difference 
form. 

B . ASSUMPTIONS 

The flow was assumed to be steady, two-dimensional and 

subsonic. For simplicity the value of specific heat (C ) was 

P 

assumed to be constant although its dependence on temperature 
and/or composition could easily be included. 

A modified Jones-Launder [8,10,11,12] two parameter tur- 
bulence model was incorporated to calculate the effective 
viscosity. It uses five empirical constants (Table I) and 
requires that two additional variables, turbulence kinetic 
energy (K) and turbuleace dissipation rate (e) , be evaluated. 
Effective viscosity was calculated using the formulas: 



The turbulent Prandtl and Schmidt numbers were taken 
equal to unity and, therefore, the turbulent Lewis nijmber was 




( 1 ) 



where 



p K“/e 



( 2 ) 



5 



c 



1 






1.43 1.92 0.09 1.0 1.3 



TABLE I. K-£ turbulence MODEL EMPIRICAL CONSTANTS 
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unity [9J. The laminar Prandtl niamber was also taken to be 
unity . 

C. GOVERNING EQUATIONS 

The conservation equations for axi- symmetrical flows 
with no tangential variations can be put into the general 
form [ 8] . 

,p.*, . i ,r^ II . i II = (3) 

convection terms diffusion terms source 

te rms 

where stands for the dependent variable (u,v,k,e,h, etc.) 

being considered (<}) = 1 for the continuity equation) , is 

the appropriate effective exchange coefficient for turbulent 
flow and is the "source term" (Table II) . The energy 

equation in terms of stagnation enthalpy has no source terms 
since the turbulent Prandtl and Schmidt numbers were chosen 
as unity and radiative transport was neglected [1,3]. Thus 
the stagnation enthalpy is given by: 

h = h + (u^ + v^)/2 + K (4) 

where for non-reacting flows: 

h = C T (5) 

P 

The calculation of temperature was made using equations 
(4) and (5). Density was calculated from the perfect gas law 





p = P/RT 


(6) 


where 


R = R/M 


(7) 
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TABLE II. GOVERNING EQUATION PARAMETERS 



D. CONSERVATION OF MASS 



On each radial line the mass flow rate was calculated 
using the local density. The error in mass flow (compared 
to the summation of "mass-in" at all upstream bovindaries) was 
used to uniformly adjust the axial velocity over the entire 
line. This process ensured that overall continuity was satis- 
fied on the line. The pressure at all downstream locations 
was then adjusted to approximately correct for the momentiam 
imbalance created by the uniform axial velocity adjustment. 

A "pressure correction" equation was then solved for each cell 
on the line. Local cell velocity (axial and radial) and pres- 
sure were then adjusted to satisfy cell-wise continuity. The 
details of this procedure are presented in reference 8. 

E. BOUNDARY CONDITIONS 
1 . Introduction 

Fixed boundaiy conditions were set at the desired or 
experimentally determined value and held constant. Specified 
gradient boundary conditions were handled by setting the 
appropriate convection/diffusion coefficient to zero in the 
finite difference equation ("breaking the link") and then 
entering the appropriate gradient through linearized "false" 
source terms [8]. The geometry and appropriate boundary con- 
ditions for the test cell are summarized in figure 1. 

2 . Inlet 

Although not a computer program limitation, "plug 
flow" was assumed at the nozzle exit and cell inlet planes. 

The (secondary) flow inlet of the TJTC was recessed approxi- 
mately 0.3 meters (figure 1) to allow a velocity profile to 
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develop over the length of the engine. Turbulence kinetic 
energy was selected to be uniform with a value which corres- 
ponded to the approximate turbulence intensity of the inlet 
flow. 

3 . Axis of Symmetry and Exit Plane 

Radial and axial gradients were set equal to zero on 
the center line and exit respectively. The radial flow velo- 
city was equated to zero. 

4 . Solid Boundaries 

All solid boundaries were considered adiabatic with 

both velocity components equal to zero ("no slip" condition) . 

For simplicity, a two part boundary layer was used. 

The border between the laminar sublayer and the turbulent 

layer was taken at y"*'=11.5 [8]. y"*” was evaluated at 

P P 

each near wall node (p) , 

where, for y > 11.5 
P - 

V = '' ■'p 

T was assumed uniform from the wall to the near wall grid 
w ^ 

point. Thus, 



y = c 
p 



1/4 ^ 1/2 

D P "^p ^/'^lam 



( 10 ) 



If yp"*” ^11.5 , the wall shear stress (t^) was calculated 
using the formula: 
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T = C p K = p C K (u/u'*’) 

w D p V// 



- ^ r n n K 1/2/ 1/4 1/2/ . 

- < Cd P \ ^p /ln(Ep5Cj5 Kp /y^^) 



( 11 ) 



where 



u* = i ln(E y/) 



( 12 ) 



Wall shear stress was evaluated for 
formula; 



+ 



< 11.5 



from the 



T 

W 



lam 



u /6 
P 



(13) 



Due to the steep gradients of properties in turbulent 
flows near solid boundaries, the source terms for K and e 
at near wall nodes were expressed in terms of the wall shear 
stress [1,8]. also provides the boundary condition for 

the u and v equations. In the following equation for tur- 
bulence dissipation rate (e) at a near wall node (p) , the 
length scale is presumed proportional to the distance from the 
wall (5) . 



e = C K ^/^/k6 = K ^/^/2.436 

P D p P 



(14) 



Equation (13) implies that the wall shear stress is calculated 

assuming a linear velocity profile when y ^ < 11.5. A near- 

P 

wall grid point, therefore, can lie within the laminar si±>layer, 
but the source terms for K and e imply that ^eff/’^lam 



11 



much greater than one [10,11]. This fact precludes from 

being significantly less than 11.5. 

F. SOLUTION PROCEDURE 

Five variables (u,v,K,e and h) were solved using equation 
(3) in finite difference form. The line by line iterative 
procedure employed upwind differencing and under relaxation to 
promote convergence [8]. Pressure (relative to a selectable 
position and magnitude within the grid) was obtained from the 
mass conservation imposed on each radial grid line and on each 
nodal control volume as discussed above. Effective viscosity, 
temperature and density were also obtained as described above. 

A more detailed explanation of this procedure can be found in 
reference 8. 
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III. TURBOJET TEST CELL NOISE CONSIDERATIONS 



A. INTRODUCTION 

Lighthill [13] and others have built a general physical 
theory of aerodynamic sound production. He has also special- 
ized that theory to the study of turbulence as a source of 
sound [14). Experiments have shown [14] that aerodynamic noise 
generated by a subsonic cold jet in the absence of any fluc- 
tuating forces between the fluid and solid boundaries exhibits 
the following characteristics: 

(1) The acoustic power output varies as a high power 
(near the eighth) of the jet velocity. 

(2) A very broad spectrum with a peak frequency near 

2 

U /2d , where U is the jet velocity and d the nozzle dia- 
meter. 

(3) Almost all the sound is radiated in directions making 
an acute angle with the jet axis. 

(a) Higher frequencies, which are apparently emitted 
mainly from the vicinity of the nozzle exhaust plane, are ra- 
diated at an angle of 45°, or slightly less, to the jet axis. 

(b) Lower frequencies, which are apparently emitted 
from the jet five to twenty diameters downstream of the nozzle, 
are radiated at much smaller angles,. 

Lighthill postulated [14] that the high frequency sound 
emanates from the heavily sheared mixing region near the nozzle 
exhaust and that the low frequency emanates from the region of 
more nearly isotropic turbulence farther downstream in the core 
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of the jet. 



While a myriad of research has been done on the generation 
and propagation of noise from a free jet, little research has 
been directed toward the problem of aerodynamic noise produced 
by confined jets. The propagation of sound associated with 
this type of flow is significantly more complicated due to the 
close proximity of several fixtures which absorb or reflect 
sound in varying degrees. In a turbojet test cell the engine 
is tested in close proximity to an augmentor txibe which is 
used for mixing and cooling the exhaust gases before they are 
vented to the atmosphere through a stack. Each of these struc- 
tures influences the sound propagation in some way. 

The numerical models discussed above can be used to predict 
velocity, pressure and turbulence intensity distributions in 
confined turbulent flows like those discussed above. It may 
be possible to use the computer codes to predict optimum place- 
ment of sound suppressors in geometrically complicated flows 
like those in a turbojet test cell. 

An attempt was made to measure the sound produced within 
the Naval Postgraduate School subscale turbojet test cell aug- 
mentor tube. These data were compared to the flow field pre- 
dictions obtained with the computer model. 

B. EXPERIMENTAL APPARATUS 
1. Test Cell 

Design and construction of the NPS subscale turbojet 
test cell are detailed in reference 15. The cell is a one- 
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eighth scale model of a NAS Alameda test cell. A TF 41 
engine was scaled to one-eighth in diameter resulting in the 
mass flow being scaled by one-sixty fourth to maintain flow 
velocities equal to those in the full scale cell. The test 
section is enclosed by hinged plexiglas sides to allow easy 
access and visual monitoring of the section during operation. 
These sides were held open during the acoustic measurements 
described in this report. The augmentor tube, equipped with 
a straight inlet, exits the cell through a removable wall. 

Its downstream end is normally attached to a deflector-plate- 
equipped vertical exhaust stack. For this experiment, however, 
the stack was removed, allowing the fluid exiting the augmentor 
tube to vent to the atmosphere. In an effort to reduce sound 
reflection inside the test cell, a one-quarter inch plywood 
divider with three-quarter inch foam padding was extended from 
the augmentor inlet to the plexiglas sides of the test cell. 

The combustor used to simulate turbojet tailpipe and nozzle 
conditions was supplied with compressed air from an Allis- 
Chalmers, twelve stage axial compressor. The tailpipe was 
fitted with either a small (2 inches diameter) or large (3.5 
inches diameter) nozzle. 

2 . Ins trument ation 

Several methods were employed in an attempt to measure 
the aerodynamic noise generated by the engine exhaust in the 
subscale test cell. The proximity of various structures in- 
cluding the augmentor tube and the stack produced an "acoustic- 
ally dirty" environment. This prevented useful data from being 
taken external to the augmentor tx±>e regardless of the pressure 
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transducer orientation to, or distance from the jet axis. The 
method chosen to obtain these measurements is illustrated in 
figure 2. 

The top of the augmentor tube incorporated a series of 
pressure taps. Several of these were modified to accept a 
one-eighth inch directional microphone. The microphone was 
inserted flush with the inside of the augmentor tiibe. Ini- 
tially it was held in a plexiglas plug with a set screw 
(figure 2) . The transducer signals were amplified, displayed 
and recorded using a Spectral Dynamics Corporation model 50330 
spectrum analyzer and a HP 7035B X-Y recorder. 

After initial testing an attempt was made to isolate those 
frequencies associated with the physical resonance of the aug- 
mentor tube and test cell from those generated directly from 
the jet mixing process. The pressure transducer was wrapped 
with styrofoam and refitted in the plexiglas holder without 
the securing set screw. To determine the degree to which the 
adjacent building was contributing to the recorded data, a 
plywood sheet was positioned between the cell and building. 

C. DATA ACQUISITION AND REDUCTION 

Sufficient data were recorded for each test condition to 
allow the mass flow rates for the nozzle, test cell and aug- 
mentor tube to be calculated. In addition, a plot of frequency 
vs sound intensity was recorded at each of the seven augmentor 
ports. On subsequent runs with the modified transducer holder 
and alternate background, data were recorded under similar 
conditions at augmentor port number four only. 
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D. COMPUTER MODEL PREDICTIONS 



At attempt was not made to predict the sound distribution 
within the augmentor tube. Rather the velocity, pressure and 
turbulence kinetic energy distributions were calculated using 
the computer program discussed above. It was desired to deter- 
mine whether these distributions could be used as a first ap- 
proximation for locating the peak noise intensities within the 
augmentor tube. 
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IV. DISCUSSION OF RESULTS 



A. INTRODUCTION . 

The puirpose of this study was to utlize a primitive vari;- 
able, finite difference computer program to analyze the flow 
within a turbojet test cell and to validate the model with data 
obtained from a subscale turbojet test cell located at the 
Naval Postgraduate School. Previous work [2,6,7] had accom- 
plished this task for a iji-aj computer model and, therefore, 
empirical data and the predictions of the model were avail- 

able for comparison. The experimental data available consisted 
of augmentor wall pressure distributions and radial velocity 
profiles along the length of the augmentor tube for low, medium, 
and high engine flow rates. 

As previously indicated, the model did a poor job of 

predicting pressure distributions in all but quiescent flow re- 
gions. Additionally, numerical convergence was difficult to ob- 
tain with that model when used to predict high velocity flows 
where compressibility effects are significant. It was antici- 
pated that a primitive variable model would help to alleviate 
these difficulties. It is desirable to have a model which can 
be used to predict the flow field for full-throttle engine con- 
ditions where the engine exhaust flow is choked. 

Adaptation of the primitive variable model to the subscale 
test cell geometry required the use of several approximations: 

(1) In modeling axi-symmetric flow, the engine was by 
necessity positioned at the axis of symmetry. In the actual 
test cell the engine was mounted closer to the deck than to 
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the overhead. It would be expected, therefore, that the 
velocity distribution in the secondary flow (the flow around 
the engine) would be somewhat different than predicted. 

(2) The subscale test cell cross section was rectangular 
while the engine and augmentor tube were cylindrical. The 
system was modeled as three concentric cylinders with cross- 
sectional areas equivalent to the physical system. The nozzle 
exit area, the test section cross-sectional area and the empir- 
ical augmentation ratios and mass flow rates were used in cal- 
culating the axial inlet velocities used in both computer models. 

(3) The actual engine incorporated a converging nozzle. 

The engine was modeled as a cylinder with a diameter equal to 
the actual nozzle exhaust diameter. 

(4) In the model the augmentor inlet and the aft test 
cell wall were taken to be flush. The actual augmentor is often 
inserted into the test section which forms a recirculation zone 
above the augmentor tube. The effects on the augmentor flow 
field introduced by this recirculation region have been shown 

to be minimal [7] . It should be noted that the model 

incorporated this recirculation zone and a flow reducing lip 
flange on the augmentor inlet. When comparisons were made bet- 
ween the predictions of the two models, the effects of this re- 
circulation zone and the augmentor inlet lip flange in the 
model were minimized by reducing their dimensions to one grid 
spacing. 

There was one additional difference between actual test 
conditions and the ntimerical model. In a test cell of speci- 
fied geometry (cross-sectional area, inlet and stack acoustic 
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and pollution treatment devices, augmentor tube length, dia- 
meter and inlet geometry, etc.) the mass flow rate of secon- 
dary air pumped by the engine exhaust and the pressure distri- 
bution throughout the cell are unique functions of the engine 
flow rate. The static pressures at the cell inlet and stack 
exhaust are both equal to the local atmospheric- pressure . 

Thus, the most desirable model capability would be the predic- 
tion of the cell flow field as a function only of engine flow 
rate. To accomplish this with a two-dimensional finite-differ- 
ence model would greatly complicate the analysis and increase 
the required computer computation time. In this investigation 
the flow field was predicted only for the test cell flow from 
the engine inlet to the exhaust of the augmentor tube. Since 
inlet boundary conditions are required for the solution of the 
elliptic equations, the mass flow rate of secondary air was 
specified at the experimentally measured valve. The predicted 
pressure profiles were then compared to the measured values. 

An alternate method could have been employed, requiring signi- 
ficantly greater computer time: iterate on the inlet mass 

flow rate until a specified augmentor exhaust pressure was 
attained. 

Three of the flow conditions selected for model valida- 
tion corresponded to conditions where empirical data were 
available. Model predictions for two additional conditions 
were made increasing the engine- to-augmen tor inlet spacing to 
one and two engine diameters. No empirical data were avail- 
able for the two latter conditions and, therefore, inlet para- 
meters were simulated using empirical data for zero engine- 
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to-augmentor entrance spacing. The test conditions are 
summarized in Table III. 

Figures 3 through 7 compare predicted axial pressure 
distributions and rddial velocity profiles obtained with the 

and u-v-p computer models. In addition, the available 
empirical data are also plotted on those figures. The velocity 
profiles from both computer models were plotted for the grid 
lines closest to the locations of the experimental data. In 
the cases where empirical data were not available, various 
representative velocity profiles were plotted for both models. 
Experimental pressure profiles were available only on the deck 
of the test cell and along the top of the augmentor tube. Thus, 
for negative positions in Figures 3 through 7 no empirical 
data were available for direct comparison to the predicted pro- 
files. The models did predict a very small pressure drop along 
the test cell deck in agreement with experiment. Predicted 
axial pressure profiles are presented for various radial posi- 
tions. These locations are given as fractions of augmentor 
tube radius (R ) . For example, the pressure distribution 

labeled R = 0.96 R indicates that the distribution is along 
an axial grid line located at a distance 96 percent of the aug- 
mentor radius from the axis of symmetry. Two model axial 

pressure profiles are depicted for each test condition. One 
profile (R = 0.38 R ) lies in the quiescent flow region bet- 
ween the engine and augmentor and was previously found by 
Walters [6] to be the only location which produced a reasonable 
estimate of the measured profiles. The second axial pressure 
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TABLE III. TURBOJET TEST CELL TEST CONDITIONS 



distribution (R = 0.13 R ) was located at less than one 

ci 

engine radius of the center line. Three u-v-p model axial 
pressure profiles are presented for each test condition. One 
distribution (R = 0.04 R ) is near the axis of symmetry. 
Another profile (R = 0.28 R ) runs along the top of the 
engine and through the turbulent mixing region. The third 
profile (R = 0.96 R ) is close to the augmentor wall. 

B. RESULTS AND DISCUSSION - PRESSURE AND VELOCITY DISTRIBUTIONS 
1 . Test Case I - Low Flow Rate/Zero Spacing 

a. Velocity Profiles 

Both models predicted virtually identical velocity 
profiles at each station along the augmentor tube (Figure 3) . 
There was close agreement between the predicted and the experi- 
mental velocity profiles at the exit of the augmentor. The 
latter result was expected since the velocity profile had be- 
come fully developed near the augmentor exit. 

b. Pressure Profiles 

(1) Model. This test condition had the 

lowest flow rates and offered the best chance for agreement 
with experiment. The initial pressure drop and the pressure 
rise were underpredicted in the outer flow region (R = 0.38 R ) 
Nearer the axis of symmetry (R = 0.13 R ) the initial pressure 
decrease was significantly overpredicted, but the magnitude 
(from minimum to maximum) and profile of the pressure rise were 
in good agreement with the experimental wall profile. The pre- 
dicted pressure curves leveled off earlier than the experimental 
data. The difficulty in obtaining good pressure profiles with 
this model are evident. Near the augmentor exit there should 
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be negligible radial pressure variations. The solution was 

converged in all dependent variables (Table V) and negligible 

pressure variation was calculated in the radial direction near 

the augmentor exit. However, the pressure profiles along the 

two axial paths (calculated by summing (9 p/8 x) . Ax from the 

inlet end) did not become equal near the augmentor exit. The 

latter resulted from the large errors in the predicted profiles 

9 D 

near the augmentor inlet (i.e., ^ in the expression 

(7/ X 

p ~ P ^ (-1^) Ax is very inaccurate near the augmentor 

X r*sr ^ dx 

x=0 

inlet) . This difficulty with the model has been previous- 

ly discussed by Walters and Netzer [6]. 

(2) u-v-p Model. Both the pressure drop and rise 
within the augmentor were underpredicted. However, the primi- 
tive variable model more accurately predicted the location of 
the minimum pressure and seemed to level off at approximately 
the same augmentor position as the experimental data. Walters 
and Netzer [6] have previously shown that the models predict 
that mixing is nearly complete at the location where the pres- 
sure profile levels off (and simultaneously, where the velocity 
distribution approaches a fully developed profile) . Experi- 
mental data confirmed this characteristic. Figure 3 indicates 
that, as anticipated, the u-v-p model can more accurately pre- 
dict the location at which turbulent mixing is complete. It 
also predicted very little pressure variation with radial 
augmentor position as is known to be true experimentally. 

2 . Test Case IIA - Medium Flow Rate/Zero Spacing 
a. Velocity Profiles 

The results are presented in Figure 4 and, again, 
both computer predictions were very similar and agreed with 
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the limited experimental data at the augmentor exit plane, 
b. Pressure Profiles 

(1) Model,. The pressure profile nearest 
the centerline became more unrealistic for this higher flow 
rate condition. The initial pressure drop was greatly exag- 
gerated; however, the magnitude of the pressure rise was again 
in good agreement with experiment. The pressure profile in 
the quiescent flow region (R = 0.38 R ) did not agree with 

Si 

the experimental curve. Both the pressure drop and rise in 
the augmentor tube were underpredicted. In addition, the 
minimtam pressure point was predicted to occur about one engine 
diameter downstream of the experimental minimum. 

(2) u-v-p Model. The u-v-p model appears to 
consistently underpredict the augmentor pressure drop and 
rise. The slope of the pressure rise, however, was in reason- 
able agreement with experiment. As discussed above, the pre- 
dicted pressure profiles were determined by specifying the 
inlet mass flow rates and uniform velocity profiles at the 
engine exhaust and cell inlet. Table IV presents the effects 
of changes in the specified flow rates on the predicted maxi- 
mum pressure rise. Also shown in Table IV are the correspond- 
ing values predicted using the one-dimensional equations for 
conservation of mass, momemtum and energy. The latter calcu- 
lations were made assuming perfect gases, coplanar engine exit 
and augmentor inlet planes, no losses, and plug flow for both 
inlets. The 1-D calculations do not correspond directly to 
the 2-D calculations (no losses and plug flow of secondary air 
at augmentor inlet) . However, the predicted magnitudes of 
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TABLE IV 

Effect of Input Variables on 
Predicted Pressure Profile 




•K 

•K 
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calculated from one-dimensional equations; co-planar inlets with 
equal static pressures and with no losses 



Ap and the effects of changing flow rates should be similar, 

max 

Table IV indicates that the 2-D and 1-D model predictions for 

the nominal case (II -A) were in good agreement. However, the 

experimental ^Pj^ax considerably larger than predicted. 

Cases II E and II F show the predicted effects of 10% changes 

in the engine and secondary flow rates, respectively. The 

predicted 1-D effects were less than the 2-D effects but in the 

correct direction and sequence. It is seen that a net 10% error 

# 

in the measured augmentation ratio (m^ from ASME orifice, 
m. from augmentor exit pitot rake, and m = m - m ) could 

w 3. u 0 

produce the differences between the experimental and 2-D model 
predictions. For example, an m^ larger by 3% and a total 
augmentor flow rate smaller by 3% (which are within the experi- 
mental accurary) would produce a decrease in augmentation ratio 
of approximately 8%. Thus, the differences between the pre- 
dicted and measured pressure profiles may not be as large as 
depicted in Figure 4. 

3 . Test Case IIB - Medium Flow Rate/One Engine 
Diameter Spacing 

a. Velocity Profiles 

Experimental data were not available for this test 
condition. The predicted profiles are presented in Figure 5 
and show that the two computer models predicted velocity pro- 
files which were in close agreement. The primitive variable 
model, however, predicted a slightly greater initial jet 
spreading at the augmentor entrance and required a slightly 
longer duct length to obtain a fully developed profile. 
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b. Pressure Profiles 

(1) Model. As for Case IIA, the centermost 

pressure distribution greatly exaggerated the initial pressure 
drop. In this case, the pressure rise also appeared to be 
much too rapid in comparison to the R = 0.38 profile 

and the u-v-p profiles. The profile in the quiescent flow 
region was in reasonable agreement with the wall profile ob- 
tained with the u-v-p model. The slopes of the pressure 
rise and the minimum pressure predicted by both models were 
nearly identical. The model, however, again predicted 

the minimum pressure to occur somewhat farther downstream than 
did the u-v-p model. 

(2) u-v-p Model. For this case, the pressure 
profile closest to the augmentor wall had a significantly 
lower minimum pressure than the other profiles. 

4 . Test Case IIC - Medium Flow Rate/Two Engine 
Diameter Spacing 

a. Velocity Profiles 

Ho empirical data were available for this test 
condition. The computer predictions are presented in Figure 6. 
As for Case IIB, the primitive variable model predicted greater 
initial jet spreading. In this case, however, the 
velocity profiles became flat and the pressure profiles leveled 
off considerably upstream of the u-v-p predicted profiles. 

b. Pressure Profiles 

(1) ip-w Model. The centermost pressure profile 
was completely unrealistic. The quiescent region profile in- 
dicated a larger pressure drop and rise and a minimum pressure 
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point farther downstream than the u-v-p model. In addition, 
the ip-u) model profile leveled off much earlier. 

(2) u-v-p Model, Again the minimum pressure was 
obtained for the profile closest to the augmentor wall and 
was located just inside the augmentor entrance. All the primi- 
tive variable curves indicate a more gentle augmentor pressure 
rise, leveling about midway down the augmentor tube. 

5 . Test Case III - High Flow Rate/Zero Spacing 
a. Velocity Profiles 

Substantially more experimental data were avail- 
able for this test condition. In the experiment the nozzle 
was operated with a pressure ratio (Pj^^^/p^) less than cri- 
tical for one-dimensional isentropic flow. However, using 
the experimental nozzle flow rate and approximating the nozzle 
flow as one-dimensional and isentropic resulted in a nozzle 
exit Mach number of approximately 0.95. This condition was 
imposed on the models. As has been observed for the previous 
test cases, the velocity profiles for both models were quite 
similar (Figure 7) and in reasonably good agreement with 
experiment. However, the predicted profiles for both models 
tended to flatten a little too rapidly. The primitive vari- 
able model again predicted slightly more initial jet spread- 
ing at the augmentor inlet and less mixing down the augmentor 
tiobe than did the ip-(o model. The experimental data more 
nearly agreed with the model at the augmentor inlet and 

with the u-v-p model downstream. 
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b. Pressure Profiles 



(1) Model. Again the center pressure 
profile was completely unrealistic. The quiescent region pres- 
sure profile underpredicted the augmentor pressure drop and 

rise. The minimum pressure position was substantially dis- 
placed down the augmentor tube and the slope of the pressure 
rise did not agree with experiment. 

(2) u-v-p Model.. All three u-v-p pressure 

profiles were quite similar. They underpredicted both the 
augmentor pressure drop and rise. The slope of the pressure 
rise was in good agreement with experimental data within the 
first half of the augmentor tube. The predicted pressure pro- 
files from both the and u-v-p models leveled off before 

the experimental curves. 

C. COMPUTATIONAL ACCURACY AND REQUIRED COMPUTER TIME 

The utility of any computer program using numerical 
methods is reflected by the amount of CPU time required and 
the ease of arriving at converged solutions. Table V compares 
the percentage change in variable magnitude on successive 
iterations and the required CPU time. 

A considerable savings in CPU time was obtained using 
the line-by-line iterative procedure of the primitive variable 
model in lieu of the point-by-point (Gauss-Seidel) method 
employed in the ijj-a) model. 

At low flow rates, the convergence was quite similar for 
both models. However, at higher flow rates the u-v-p model 
had better convergence in less time. 
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TABLE V. COMPARISONS OF CPU TIME AND (A<l)/(t)) • 100 FOR LAST 
CALCULATION 



D. COMPUTER RELATED PROBLEMS 



Both models required that the proper relaxation para- 
meters be selected in order to obtain convergence. The lack 
of any procedure for selecting the proper relaxation values 
makes this process quite time consuming. Previous research 
using the model facilitated the selection of these 

parameters for that model. It was found that the u-v-p model 
was quite sensitive to the calculated "pressure corrections". 
Obtaining the correct underrelaxation value for pressure 
proved to be the key in arriving at a converged solution for 
the primitive variable model. 

The line-by-line iterative procedure used in the u-v-p 
model was, as expected, quite good in propagating distur- 
bances downstream when sweeping from left to right. A down- 
stream disturbance is propagated upstream by successive sweeps 
through the entire field. This fact, at least for a geometry 
incorporating a sudden contraction, made the convergence de- 
pendent on the number of traverses on each radial line. An 
excessive number of traverses would cause divergence. The 
number of traverses on each line was controlled in two ways. 
After each traverse, residual factors were calculated for 
each variable and the largest residual factor was compared to 
a pre-set value. If the largest residual was less than the 
pre-set value, the program advanced to the next radial line. 
The program would also advance when a pre-set maximum number 
of traverses had been completed on any radial line. To aid 
convergence, only a few traverses on each line were allowed 
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until several sweeps through the field had allowed the pre- 
sence of the contraction wall to be "felt” upstream. The 
number of traverses on a line was then increased. It was 
additionally found that when working with coaxial flows with 
radically different inlet velocities, the normalizing factors 
(which were based on average inlet conditions and used to 
calculate the residual values on each line) resulted in ex- 
cessive traverses being made in regions of high flow velocity. 
Repeated calculations on radial lines which had already con- 
verged often caused divergence. Adjusting the normalizing fac- 
tors downstream of the engine exit alleviated this problem. 

The primitive variable model demonstrated some conver- 
gence difficulty in the recirculation region adjacent to the 
sudden contraction. This problem could have resulted from 
the relatively large normalizing factors used in this local 
region of low velocity or, as suggested by Launder and 
Spalding [12], it could have possibly been due to the in- 
adequacy of the empirical constants (Table I) in the K-e 
model for the flow in this quiescent zone. 

As with any finite difference numerical solution, grid 
spacing was found to be critical. To aid convergence, the 
grid spacing was decreased in all regions where large pro- 
perty gradients were expected. Gosman, et al [1] recommended 
that for the model, successive spacing should not in- 

crease by more than about a factor of 1.5. This restriction 
was also employed for the primitive variable model. A 30 by 
30 grid system was utilized for the primitive variable model. 
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whereas a 43 by 40 grid was required for the model. 

E. SUMMARY OF RESULTS - PRESSURE AND VELOCITY DISTRIBUTIONS 

In most cases the predicted velocity profiles for both 
the ip-u) and u-v-p computer models agreed with each other and 
with the available experimental data. The u-v-p model predic- 
ted more initial jet spreading at the augmentor inlet, espe- 
cially in those cases where the engine exit plane was not flush 
with augmentor inlet. Both models predicted that a flat velo- 
city profile was obtained where the pressure profile leveled 
off. For the one case in which experimental data were avail- 
able over the entire legnth of the augmentor tube, the pre- 
dicted velocity profiles seemed to flatten slightly faster 
than the experimental data. 

In general, the ij;-w model demonstrated poor pressure 
prediction capability. For low flow rates the pressure rise 
on the wall from minimum to maximvim was accurately predicted 
for a pressure profile calculated near the axis of symmetry. 

For higher flow rates the centermost predictions became un- 
realistic although converged solutions were obtained for all 
primary variables. For the predicted pressure profiles along 
axial lines that were in the quiescent flow region, the ip-ui 
model characteristically unpredicted the pressure drop and 
rise along the augmentor wall. The minimum predicted pres- 
sure location was typically somewhat downstream of its experi- 
mentally determined position. 

The primitive variable model predicted little radial pres- 
sure variations except near the engine exit plane. It seemed 
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to accurately locate the minimum pressure position and to pre- 
dict the rapid pressure drop at the augmentor entrance. It 
consistently underpredicted the pressure drop and rise in the 
augmentor but tended to accurately predict the slopes of the 
rising pressure profiles. However, as discussed above, reason- 
ably expected experimental errors could account for these dif- 
ferences. Application of the K-e turbulence model with 
fixed parameters (C , C,, C_, a , a„) to the test cell flow 
conditions may also be a reason for the lack of quantitative 
accuracy in the predicted pressure profiles. Further experi- 
mental data (for example, turbulence intensity measurements) 
are needed to further validate the model. 

The primitive variable model appears to reasonably pre- 
dict the flow field up to sonic engine exhaust conditions. 

These predictions include realistic pressure distributions 
and require substantially less computer time and fewer grid 
lines than the ip-cj model. 

F. RESULTS AND DISCUSSION - NOISE CONSIDERATIONS 

Computer model predictions for the velocity, turbulence 
kinetic energy and pressure distributions are presented in 
Figures 8, 9 and 10. The predictions were for the two 
inches diameter nozzle at high mass flow rate (Run 1, Table 
VI) . The four curves on each figure are for different radial 
positions; jet centerline, just inside the jet exhaust radius; 
just outside the jet exhaust radius, and near the augmentor 
wall. 
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TABLE VI. SUMMARY OF TEST DATA FOR NOISE EXPERIMENTS 



Noise intensities for selected frequency ranges as a 
function of augmentor txibe position are displayed in Figures 
11, 12 and 13. 

Figures 14 through 16 and Table VII compare the basic 
test results to data obtained at one transducer port after 
the microphone was insulated and the "building wall" location 
was changed. 

Figure 11 presents the data for the small nozzle with a 
' high svibsonic exhaust velocity (M^ a; 0.75). The highest 
frequencies were greater than peak values expected from ori- 
fice data (V/2d) but were apparently generated near the 
nozzle exit as expected. If the high frequencies generated 
near the nozzle exit were radiated at 45° to the axis, the 
corresponding peak intensities would occur between two and 
five jet diameters downstream. The lowest frequencies peaked 
in intensity at approximately 20 jet diameters downstream in 
general agreement with data from free jets emitted from orifice 
plates [14]. The intensity increased as the frequency de- 
creased. 

Figure 12 presents the data for the small nozzle with a 
V very low exit velocity (M^ ~ 0.25) . The data had similar 
characteristics to that for the high exit velocity: (1) the 

intensity increased with decreasing frequency, (2) the peak 
frequencies were greater than V/2d and (3) the lowest fre- 
quencies peaked in intensity at approximately 20 jet diameters 
downstream. In this case, however, the highest frequencies 
peaked in intensity further downstream. 
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TABLE VII. TEST DATA FOR MODIFIED NOISE EXPERIMENTS 



Figure 13 presents data for the larger nozzle with low 
exit velocity. The data were similar to Figure 12 except 
that the intensity did not vary as consistently with frequency 
range. The peak intensity did occur for the lowest frequency 
range and occurred at approximately 10 jet diameters down- 
stream. This peak intensity was at the same augmentor posi- 
tion as for the smaller nozzle at low exit velocity (Figure 12) . 

The possibility existed that some of the sound recorded 
was either being reflected from or was due to the physical 
resonance of surrounding structures. However, the data pre- 
sented in Figures 14 and through 16 indicate that these effects 
were not large. 

The computer predictions made for the distributions of 
axial velocity, turbulence kinetic energy and pressure (Figures 
8, 9 and 10) were for the test conditions/results presented in 
Figure 11. (RUN No. 1, Table VI.) 

The peak-intensities for the lower frequency sound noted 
above occurred approximately at the predicted location within 
the augmentor tube at which, (1) turbulent mixing became essen- 
tially complete, (2) pressure reached a maximum, and (3) the 
axial velocity profile became fairly constant. The decay in 
the predicted turbulence kinetic energy (Figure 9) followed 
closely the rise in noise amplitude for the two lower fre- 
quency ranges (Figure 11) . 

The data were much too limited to draw any general con- 
clusions. However, the results indicate that the maximum 
sound intensities occur within the augmentor tube in the 
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region where jet mixing is nearing completion. If this ob- 
servation is correct, it would enable the computer model to 
be used for full scale test cells to predict the best place- 
ment of sound suppressors/jet breakup devices for varying 
test cell geometries and engine operating conditions. 
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V. CONCLUSIONS 



The primitive variable computer model for the turbojet 
test cell appears to reasonably predict the flow field for 
engine exhaust conditions up to the sonic condition. The 
conclusion is based upon subscale test cell measurements of 
velocity and pressure distributions. To further validate the 
model# comparisons should be made with full scale test cell 
data. In addition, turbulence intensity distributions mea- 
sured within the subscale test cell would provide valuable 
new data for additional model validation. 
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(DIMENSIONS IN INCHES) 




CO 

o 

0. 

cc 

o 

h- 

2 

UJ 

2 

CO 

Z) 

< 



46 



FIGURE 3. AUGMENTOR PRESSURE AND VELOCITY PROFILES FOR STRAIGHT 

PIPE INLET AND ZERO ENGINE -AUGMENTOR SPACING, = 0-45 kg/sec 
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FIGURE 4. AUGMENTOR PRESSURE AND VELOCITY PROFILES FOR STRAIGHT 

PIPE INLET AND ZERO ENGINE-AUGMENTOR SPACING, A = 0.69 kg/sec 
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FIGURE 6. AUGMENTOR PRESSURE AND VELOCITY PROFILES FOR STRAIGHT 
PIPE INLET AND TWO DIAMETER ENGINE-AUGMENTOR SPACING, 
m =0.69 kg/sec 
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FIGURE 7. AUGMENTOR PRESSURE AND VELOCITY PROFILES FOR STRAIGHT 

PIPE INLET AND ZERO ENGINE-AUGMENTOR SPACING, A = 0.99 kg/sec 
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FIGURE 8. PREDICTED VELOCITY PROFILES 
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FIGURE 11. NOISE INTENSITY VS. AUGMENTOR TUBE POSITION, RUN NO 
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FIGURE 12. NOISE INTENSITY VS. AUGMENTOR TUBE POSITION, RUN NO. 
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FIGURE 13. NOISE INTENSITY VS. AUGMENTOR TUBE POSITION, RUN NO. 
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FIGURE 14. NOISE INTENSITY VS. FREQUENCY, RUNS 1 AND la 



TEST CONDITIONS 
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FIGURE 15. NOISE INTENSITY VS. FREQUENCY, RUNS 2 AND 2a 
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FIGURE 16. NOISE INTENSITY VS. FREQUENCY, RUNS 3 AND 3a 
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